Spin Bose Glass Phase in Bilayer Quantum Hall Systems at u = 2 
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We develop an effective spin theory to describe magnetic properties of the v — 2 Quantum Hall 
bilayer systems. In the absence of disorder this theory gives quantitative agreement with the re- 
sults of microscopic Hartree-Fock calculations, and for finite disorder it predicts the existence of a 
novel spin Bose glass phase. The Bose glass is characterized by the presence of domains of canted 
antiferromagnetic phase with zero average antiferromagnetic order and short range mean antiferro- 
magnetic correlations. It has infinite antiferromagnetic transverse susceptibility, finite longitudinal 
spin susceptibility and specific heat linear in temperature. Transition from the canted antiferromag- 
net phase to the spin Bose glass phase is characterized by a universal value of the longitudinal spin 
conductance. 
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PACS numbers: 73.40. Hm, 73.20.Dx, 75.30. Kz 

Recently a canted antiferromagnetic phase has been 
predicted in bilayer quantum Hall (QH) systems at a 
total filling factor v — 2 on the basis of microscopic 
Hartree-Fock calculations and a long wavelength quan- 
tum 0(3) nonlinear sigma model In this letter we 
construct an alternative effective spin theory that can 
describe the richness of the phase diagram of a bilayer 
u = 2 quantum Hall system. Our effective spin the- 
ory treats the interlayer tunneling nonperturbatively, in 
contrast to the 0(3) nonlinear sigma model which in- 
cludes tunneling perturbatively through an antiferromag- 
netic exchange. It gives excellent agreement with the re- 
sults of microscopic Hartree-Fock calculations in and 
extends the earlier effective field theory by allowing to 
study quantitatively the effect of a finite gate- voltage be- 
tween the layers and calculate intersubband excitation 
energies. Our theory can easily incorporate the effects of 
disorder and we predict that for any non-zero disorder 
there is a new v = 2 spin Bose glass quantum Hall phase 
which may be visualized as domains of canted antiferro- 
magnetic phase surrounded by domains of fully polarized 
ferromagnetic or spin singlet phases. In this system the 
Bose glass phase we predict is quite novel, and we elabo- 
rate in this Letter on the origin and the properties of this 
new QH glass phase. Related disorder induced spin phase 
has been discussed in a different setting in reference |Q. 

In the absence of interlayer interaction each layer of 
the V — 2 bilayer system would be in a fully spin polar- 
ized ferromagnetic v = \ incompressible QH state with 
spins in both layers pointing in the direction of the ap- 
plied magnetic field (FPF state). Tunneling between the 
layers favors the formation of spin singlet states from 
the pairs of electrons in the opposite layers and energet- 
ically stabilizes the spin singlet (SS) state. In [|| it was 
observed that the competition between the two tenden- 
cies may lead to a third intermediate phase: canted an- 
tiferromagnetic state, where spins in the two layers have 
the same component along the applied field but opposite 
components in the perpendicular 2D plane (CAF state). 



We now introduce a simple lattice model which we use 
to describe the physics of the bilayer v = 2 QH system. 
We consider a bilayer lattice model shown in figure ||. 
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FIG. 1. Effective bilayer lattice model for v = 2 double 
layer QH system 

Sites in each layer may be thought of as labeling dif- 
ferent intra-Landau-level states. Electrons may tunnel 
from one layer to another conserving the in-plane site in- 
dex (i.e. between the states with the same intra-Landau- 
level index). There is a ferromagnetic interaction be- 
tween nearest-neighbor sites within individual layers and 
a Zeeman interaction with the applied magnetic field. We 
also account for the charging energy, i.e. the energy cost 
of creating charge imbalance between the layers through 
the term Tic below. The Hamiltonian of the system may 
be written as 

7i — Ht + Tic + T^z -\-'Hf 
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where T or _B is the isospin index that labels electrons 
in the top and bottom layers respectively, i is the in- 
plane site ( intra-Landau-level ) index, and a is the spin 
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index. S'^- = dp^^a^pCnp and nxi = Y.a'^Tia'^Tia are 
spin and charge operators for layer T, with analogous 
definitions for layer B. Parameters J and ec of this model 
may be easily estimated as p'g = e^/(16V27reZ) , where 
I = \fhcJeB is magnetic length and tc = ^\f\^ ~ 
e~'^^ ^^'^ Er f[d/ {1^/2)]) , where d is the distance between 
the layers and Erf is the error function Q . 

Let us consider an individual nmg, i.e. two sites with 
the same in-plane site index on the opposite layers. Each 
rung must be populated by two electrons, therefore we 
have six possible states for each rung. The six available 
states are conveniently classified into three states that 
are spin triplets 
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and three states that are spin singlets 
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Operators t and r satisfy bosonic commutation relations 
and constraint r^Ta+tl^ta = 1 projects into the phys- 
ical Hilbert space. 

In Hamiltonian (|l]) all the terms except Hf act within 
a single rung. It is therefore natural as a first step to 
diagonahze H' = Ht +'Hc+ Hz on one rung. The latter 
task is simplified by the observation that Ht and He 
act in the subspace of r states, whereas Hz operates in 
the t subspace. A simple calculation gives for the lowest 
energy eigenstates of H': 

• state |t+) with energy Et ^ —Hz 



• state 
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with energy £:„ = ^ - V'^sAS + ec/4- 
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State is a spin-singlet state whose energy is lowered 
by interlayer tunneling and is a spin triplet state 
favored by Zeeman interaction. Competition between the 
two states is a competition between the SS state and the 
FPF state. In the absence of the in-plane ferromagnetic 
interaction we would have level crossing at Ey — Et with 
a first order phase transition between SS and FPF phases. 
However as we show below Hp acts as an interaction that 
connects the two states and gives rise to an intermediate 
state that is a superposition of the \v-[.) and states 
and corresponds to the CAF phase. 

We rewrite Hamiltonian (|l|) keeping only the lowest 
energy states and 



H = EtY^ t\^U+ + E, v\^v,+ 

i i 

j{cos0 + sin 61)2 (tl^v,+ vj^tj+ + t]^Vj+vj^ti+^ 



and the hard core constraint is implied 
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The mean field analysis of Hamiltonian may be 
done by 
considering states with simultaneously condensed v and 
t bosons. They correspond to the variational wavefunc- 
tions of the form |*) = exp{a J2^ vl+ + PE^ 4+}|0> §■ 
The energy of state \'^) is given by 



Eo = Ey\a\^ + EtW - Jicosd + sirL9y\a 
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and state I^E") obeys constraint (||) on the average pro- 
vided that 



a- 



1 



(7) 



Values of a and /3 that minimize (g) under the condition 
(0) are given by 



|a| = i 
lal =0 
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where 



J(cos 9 + sin 6)^ - Et + E^ 
2[J(cos6'-Fsin6')2 - J] 
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The first and the last cases obviously correspond to the 
FPF and SS phases respectively But there is also a 
nontrivial new phase that appears in our analysis when 
both \a\ and \(3\ are finite. It is easy to verify that this 
state corresponds precisely to the canted antiferromag- 
netic phase discussed in 0| with direction of the Neel 
ordering given by the phase between the and v+ con- 
densates tan'^Ny/Nx — Arg{a* 0) ( Neel order parame- 
ter is defined as N = Sn — X)i ^Bi )■ On figure || we 
show the phase diagram obtained from equation (H) . We 
can also use our bosonic model to calculate the phase di- 
agram in the presence of an interlayer charge imbalance 
and these results will be reported elsewhere 

The lowest energy interband transition in the SS phase 
will correspond to destroying a f + and creating a bo- 
son. The energy for such transition is w_ = -~J{cos6 + 
sind)'^ + Et — Ey and vanishes at the SS/CAF transition 
as may be seen from equations and (^. Analogously 
in the FPF state the lowest energy interband transition 
will correspond to destroying and creating a boson 
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FIG. 2. Phase diagram for disorder-free system. All ener- 
gies are in units of J, Ec = 5.0J. 

Let us now give a simple physical picture that will illus- 
trate formal calculations presented above. We consider 
an SS state which has singlet w^-bosons on all rungs and 
imagine creating i+-triplet on one of the rungs. Creating 
a localized triplet requires energy Et — and this en- 
ergy is unaffected by the ferromagnetic interaction since 
parallel and anti-parallel contributions cancel for triplet 
interacting with neighboring singlets. However Tip also 
gives rise to a process in which one of the spins of the 
triplet pair and one spin from the neighboring singlet 
pair are flipped simultaneously This process is shown on 
figure |3| and may be interpreted as hopping of the triplet 
boson to the nearest-neighbor site. Therefore creating a 
propagating triplet boson at wavevector k will give it an 
additional kinetic energy J{coskx + cosky) due to Tip- 
This allows us to have a situation when Et ~ E^ > but 
Et — Ey — 2J < 0, i.e. when it is energetically unfovarable 
to create localized triplets but it is already favorable to 
create them at fc = 0, i.e. to have a condensate of t+ 
bosons. This effect is the origin of the CAF state and 
allows us to understand this phase as a coherent super- 
position of condensed t+ and bosons. 
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FIG. 3. Triplet "hopping" process due to ferromagnetic 
in-pane interaction 

In a real system there is always disorder. It may be 
due to fluctuations in the distance between the wells or 
the presence of impurities. Such disorder can be easily 
included in our effective bosonic theory, but would be 
difficult if not impossible, to include in the Hartree-Fock 
theory of references 0| . For our effective spin model the 
major effects of disorder will be randomness in the value 
of tunneling Ag^^s and the appearance of a random lo- 
cal gate voltage, in both cases leading to random local 
fluctuations in the energy of the w+ boson. Then if we 



are close to the CAF-SS transition we may have a sit- 
uation induced by disorder where Et — E™"-^ — 2 J < 
and Et — i?™™ — 2J > 0. So for some regions creating 
non-local t+ triplets will lower the energy of the system 
and for some regions it will lead to an energy increase. 
In this case the system breaks into domains, with each 
domain being locally a CAF phase or a SS phase (region 
III on figure Each CAF domain may be thought of as 
being in a quantum disordered state with an undefined 
direction of the Neel order but finite z-magnetization . 
Close to the CAF-FPF transition line in the disorder-free 
system we may have a disorder induced situation where 
we have CAF domains in the background of domains of 
the FPF phase (region I on figure ^). Finally we can also 
have the phase where we have domains of all three kinds 
( region II on figure || ) . In figure || we show the result- 
ing phase diagram for the same values of parameters as 
in figure || but assuming that As^s may randomly vary 
by 10% around its average value. Such a variation in 
AsAS is physically quite reasonable even in high quality 
2D systems since As as depends exponentially on layer 
thickness. There are no phase transitions between re- 
gions I, II and III on the phase diagram in figure |^ but 
only smooth crossovers. The true quantum phase transi- 
tions occur between FPF and I, SS and III and between 
CAF and one of the I, II or III regions. 
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FIG. 4. Phase diagram for a disordered system. All en- 
ergies are in units of J, ec = 5 J. Variation in Asas was 
assumed to be 10 %. Region I corresponds to domains of 
CAF and FPF phases. Region III to domains of SS and CAF 
and region II to domain of all three kinds. There are no phase 
transitions among regions I, II, III, they all correspond to the 
spin Bose glass phase. 

The nature of these phase transitions is also easy to 
understand. The SS phase is an insulating phase of zero 
density of t+ bosons, FPF state is an insulating phase 
with density n = 1 Q and CAF is a superfluid phase. 
Randomness that we consider acts as a randomness in 
the chemical potential of these t_|_ bosons, so our prob- 
lem is equivalent to the problem of bosons in a random 
potential, the so-called dirty boson problem, considered 
in references [^,0. We immediately recognize I, II and 
III as a single Bose glass (SBC) phase of the singlet and 
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triplet bosons. This observation allows us to draw sev- 
eral important conclusions about the properties of this 
SBG phase. In the SS state the (S'^(-w)5'^(w)) corre- 
lation function is zero and in the CAF phase it has a 
(5-function peak at zero frequency due to the Goldstone 
mode of the spontaneous breaking of the U{1) symmetry 
of spin rotations around the z-axis. In the SBG phase 
this correlation function will be finite at small frequen- 
cies, which implies finite longitudinal spin susceptibility 
and is the analog of finite compressibility of the usual 
charge Bose glass. Our new SBG phase does not have 
antiferromagnetic long range order, i.e. {N^i^y)) — 0. All 
the i+ and bosons are localized in this phase, therefore 
it will have only short range mean antiferromagnetic cor- 
relations. But analogous to the infinite superfiuid suscep- 
tibility of charge Bose glass our SBG phase will have an 
infinite transverse antiferromagnetic susceptibility. An- 
other important feature of the Bose glass phase is a finite 
density of low energy excitations [^,0 . This implies that 
our SBG phase will have a specific heat linear in tem- 
perature which provides another way to experimentally 
distinguishing it from the CAF phase whose specific heat 
goes as or FPF and SS phases that have exponen- 
tially small specific heat at low temperatures. The ex- 
istence of the SBG phase separating SS, FPF and CAF 
phases also has important consequences in that it changes 
the critical exponents for the corresponding phase tran- 
sitions from the one obtained in ^ for the disorder-free 
system. The new critical exponents will be those of the 
superconductor-insulator transition in dirty boson sys- 
tem studied in [^j^ . We would also like to point out that 
the SBG system that we suggested may be a better ex- 
perimental realization of a 2d superconductor-insulator 
transition in a boson system than conventionally used 2d 
superconducting films § in that it is free of long range 
forces and allows one to vary the density of bosons by 
varying H^- In addition our predicted QH Bose glass 
phase transition does not have the complication arising 
from parallel fermionic excitations which may play a role 
in the superconducting films We therefore expect, 
based on the arguments given in j6|Q|, that transition 
from the CAF phase to the SBG phase will be character- 
ized by a truly universal longitudinal spin conductance, 
which in principle can be measured by measuring the spin 
susceptibility and the spin diffusion coefficient. 

Before concluding we remark on the feasible exper- 
imental observability of our proposed QH Bose glass 
phase. First we remark that the basic v = 2 QH phase 
transition and the associated softening of the relevant 
spin density excitations has been verified experimentally 
via inelastic light scattering spectroscopy. Since in- 
terlayer tunneling fluctuations are invariably present in 
real systems, it is in fact quite possible that the exper- 
iments in have already observed a transition to the 
Bose glass phase as in our figure ^ Some evidence sup- 
porting this possibility comes from the fact that softening 



of the spin density excitations observed in |^ did not 
lead to the appearance of a sharp dispersing Goldstone 
mode expected in the CAF phase but only to some broad 
zero energy spectral weight consistent with the Bose glass 
phase. Future experiments in samples with deliberately 
controlled disorder should be carried out to conclusively 
verify our prediction of a QH disordered Bose glass phase. 

In conclusion we predict a new 2D Bose glass phase 
in a = 2 QH bilayer system by introducing an effective 
spin theory. This phase has the usual properties of a Bose 
glass phase |^J^ including a universal spin conductance 
at the transition. While we have specifically considered 
the V = 2 integer QH situation, our arguments should go 
through for all v = 2/ (odd integer) fractional QH states 
also, following the reasoning of Q, and for the fractional 
filling there should be an exotic fractional quantum 2D 
Bose glass in bilayer systems ||TT| . 
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